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Abstract 

We consider triplet Higgs model with A4 symmetry to generate the neutrino mass 
matrix. The tribimaximal form of the neutrino mixing matrix can be naturally ob- 
tained. Imposing the neutrino oscillation data, we show that 1) both normal and 
inverted mass hierarchy are allowed, 2) there is a lower bound on the lightest neutrino 
mass and the effective mass for neutrinoless double beta decay, 3) the non-vanishing 
#13 can be accommodated by considering small perturbation, 4) # a tm should be very 
close to 7r/4 even after perturbation. 
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1 Introduction 



During last decade we have seen the firm evidences that the neutrinos have masses and 
mixings which may indicate the new physics (NP) beyond the Standard Model (SM). 
The solar, atmospheric and reactor neutrino experiments have measured the mass squared 
differences and mixing angles [1] as shown in Table 1. From the data in Table 1 only, 



Parameter 


Best fit 


3-0" c.l. range 


Am 2 (10~ 5 eV 2 ) 

AmL n (l(T 3 eV 2 ) 
sm 2 (# ) 

sin 2 (6» atm ) 
sin 2 (# 13 ) 


(7.65183) 
(2.40t°;} 2 ) 

0.3041™ 22 . 

U - OU -0.06 

n m+ - 016 

u - Ui -o.on 


7.05 - 8.34 
2.07-2.75 
0.25-0.37 
0.36 - 0.67 
< 0.056 



Table 1: The neutrino oscillation data. 



we do not know the mass hierarchy, the absolute neutrino masses, or whether neutrinos 
are Majorana or Dirac particles. The future experiments like the neutrinoless double beta 
decay or Tritium decay may give further information on the nature of neutrinos. Therefore, 
it is important to get a model which predicts the mass matrix of neutrinos. 

The current experimental data for the neutrino oscillations suggest that the mixing ma- 
trix, the Pontecorvo-Maki-Nakagawa-Sakata (MNS) matrix, is tribimaximal at the zeroth 
order [2]: 

( 2/3 1/3 ^ 
(\U lv \ 2 ) = 1/6 1/3 1/2 • (1) 
V 1/6 1/3 1/2 j 

One of natural flavor symmetries which give the tribimaximal MNS matrix is the A 4 sym- 
metry [3, 4]. In addition, it was shown that the effective mass parameter relevant to the 
neutrinoless double beta decay can be predicted in a model with A 4 flavor symmetry [4]. 
The tribimaximal mixing can also be naturally realized in the triplet Higgs model. 
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The triplet Higgs model, a TeV scale NP model, can generate Majorana neutrino masses. 
With symmetry the Higgs triplet model naturally gives tribimaximal mixings to neutri- 
nos. The predicted doubly charged Higgs (if ±=t ), if light, can be found in the CERN Large 
Hadron Collider (LHC) experiment [5]. 

In this Letter we study a triplet Higgs model with A4 flavor symmetry which predicts 
not only the neutrinoless double beta decay but also all the neutrino masses. The Letter is 
organized as follows: In Section 2 we introduce our model. The numerical analysis for the 
neutrino masses and mixings and the predictions for the neutrinoless double beta decay 
are done in Section 3. We conclude in Section 4. 



2 The Higgs triplet model with symmetry 

The conventional method to generate the neutrino masses is to introduce the heavy Ma- 
jorana right-handed neutrinos (type-I seesaw). However, it is difficult to test this case 
experimentally. An alternative way to give Majorana masses to neutrinos is to introduce 
SU(2) L Higgs triplet A with U(l)y charge 1 in the SM (type-II seesaw) [8]. Via Yukawa 
interaction the Higgs triplet model can provide Majorana masses to neutrinos if the neutral 
component A gets very small vev [7]. The Yukawa interaction is given by 

C Y = ^LjCiT^Lj + H.c, (2) 

where hij is a complex symmetric coupling matrix, Lj = (i/j, is a left-handed lepton 
doublet, C is the Dirac charge conjugation operator, and r 2 is a Pauli matrix. The Higgs 
triplet can be decomposed as follows: 

a = / A+/V2 A++ 
^ A -A+/y/2 
Then the neutrino mass matrix is written in terms of the vev as 



(3) 



{M v )ij = {A). (4) 

There are many attempts to obtain the form of neutrino mass matrix suggested by 
Table 1 in the framework of A 4 symmetry in the literature [3, 4]. The group A 4 has three 
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1" 


3 
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-1/2 
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Table 2: The assignments of A 4 and the SU(2) L x U(l)y representations in our model. 

singlet representations, 1, 1', 1", and a triplet representation, 3. Their tensor products are 
decomposed as 

3©3 = 3 s ©3 a ©l©l'© 1", 1' © 1" = 1. (5) 

For two vectors (xi, x 2 , X3) and (1/1,1/2,1/3) transforming as 3, the first rule in the above 
equation states that 

(3©3) 3s = (x 2 y 3 + x 3 y 2 ,x 3 y 1 +x 1 y 3 ,x 1 y 2 + x 2 yi), 

(3 © 3) 3a = (x 2 y 3 - x 3 y 2 , x 3 y x - x ± y 3 , x x y 2 - x 2 y x ), 

(3 ©3)i = x x yx + x 2 y 2 + x 3 y 3 , 

(3 ©3)i/ = xiyi + ujx 2 y 2 + uj 2 x 3 y 3 , 

(3©3)i» = xiyi + u 2 x 2 y 2 + ux 3 y 3 . (6) 

We assign the three-dimensional representation 3 of A 4 to the doublet Higgs <fi and the 
triplet Higgs A as in Table 2. To implement flavor symmetric Lagrangian and get a 
neutrino mass matrix consistent with the experiments, we need to introduce additional 
Higgs multiplets. Here we introduce three more triplet Higgs x% (i — 1,2,3) like [4], which 
are assigned 1, 1', 1" of A4. We assign 1, 1', 1" to the right-handed leptons and 3 to the 
left-handed leptons. 

Then not all the couplings in (2) but only the following A 4 -symmetric Yukawa interac- 
tions are allowed: 

C = -A e e i? (0 t L) 1 -A / e ? 7 i? (0 t L) r -A / e / ^(0 t L) 1 » 

+ ^\ a L t Cit 2 AL + X -\ x {L T Cir 2 XxL) x + 1 -\ 2 {L T Cit 2 x 2 L) v , + ^X 3 (L T CiT 2X3 L){7,) 
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where the subscripts 1,1', 1" represent the transformation rules of the ((f)* L) pair in the 
first line and of the (L T L) pair in the second line. The resulting lepton and neutrino mass 
matrices have the following form 



Mi 



M,, 



\ e V\ X e V 2 X e V3 
\' e Vi u\' e V2 UJ 2 \' e V3 
\^ \"V! U 2 \"V 2 U}\"V 3 J 

t a+b+c f 

f a + u 2 b + ojc 
ye d 



(8) 



e 
d 

a + ub + u 2 c 



\ 



(9) 



where Vi = (0, 



1, 2, 3) and 



a = Ai(x?>, b = \ 2 { X l), c = \ 3 ( X ° 3 ), 
d = X A (A° 1 ), e = A A (A°), / = Aa(A 



(10) 



We assume Vi = v 2 = v 3 , and d = e = f, which can be guaranteed by a residual symmetry 
like Z 3 . For simplicity we impose additional assumption: b = c which naturally gives the 
maximal atmospheric neutrino mixing, although it is not required by any symmetry. It is 
straightforward to extend our analysis to the case b ^ c. We will consider the effect of 
small perturbation of c — b in Section 3.2. 

Then the lepton mass matrix, M h can be diagonalized by rotating the left-handed lepton 
by the unitary matrix 

fi i i \ 

1 uj u 2 
1 u 2 



V 



The neutrino mass matrix is diagonalized by the transformation 



U T v M v U v = Mf^, 



(12) 



where U u is a unitary matrix which can be decomposed into three successive rotation 
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matrices: U v = V23V13V12 [1]. The form of the neutrino mass matrix when b = c is given by 

/ r. _L Oh A A \ 



M,, 



\ 



a + 2b 


d 


d 


d 


a — b 


d 


d 


d 


a — b 



(13) 



The resulting MNS matrix is the product of U(u) and U u : 

Umns = U(uj) U v . 



(14) 



Alternatively we can work in the basis where the mass matrix of the charged lepton is 
diagonal. In this basis the neutrino mass matrix is in the form: 



m: 



cff 



U(uYM u U{u) ] 



( a + 2d b b ^ 
b b a — d 

b a — d b 



(15) 



such that 



UT s MfU 



MNS lvl u U MNS 



M, diag . 



(16) 



Since the neutrino mass matrix is symmetric under the exchange of -h- u t , we get a 
bimaximal 2 — 3 mixing and a vanishing 1 — 3 mixing 



23 = 7r/4, e 



13 



0. 



;i7) 



The (15) does not automatically generate the desired mixing angle #12 for the tribimaximal 
mixing matrix. The form of (1), however, can be achieved in a wide region of parameter 
space. This can be seen from the fact that the Umns obtained in (16) can be decomposed 
in general as 



C/mns = diag(e^,e^,e^) f^ NS = P(/« NS , 



(18) 



where ?7^ NS is the physically measurable matrix which includes the Majorana phases in 
the standard form [1], 

/ C12C13 s 12 c 13 s n e~ lS \ 



77 st 



-S12C23 - Ci2S2 3 Sl 3 e C12C23 - Si2S 23 Sl3e S23C13 

V S 12 S 2 3 - Ci2C 23 Sl3e 1 ' 5 -C12S23 - Sl2C23Sl 3 e l<5 C23C13 J 
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x diag(e iai , e ia \ 1).(19) 



The phase matrix, P(/3), in (18) can be absorbed into the right-handed charged lepton 
sector, and therefore, is unphysical. However, since the is complex matrix in general, 
the phase matrix P(/3) is generally allowed, and it makes the allowed parameter space much 
more wider than the real matrix case. The condition for the tribimaximal mixing including 
the unphysical /3's can be written as 

( e tfi _ e ^ 2 )[a(e ift + e i/?2 ) + fee^ 2 )] 

where the phases f3±, f3 2 can take arbitrary values. The d = is a special case to give the 
tribimaximal mixing matrix. 

We note that if we assume (A°) = (A3) = 0, the tribimaximal form can be obtained 
without further conditions among a,b,d. This case was considered in [6]. 

Our model is similar to the one considered in [4]. But their neutrino mass matrix is 
different from ours and can be obtained by exchanging b ■<->■ d. Phenomenologically the two 
are significantly different in that, for example, d = gives the tribimaximal mixing matrix 
in our case, while the model in [4] cannot. 



3 The numerical analysis and model predictions 



3.1 The case b = c 

As mentioned in the previous section, we get 6*23 = 7r/4, 9 13 — in this case. More explicitly, 
with U MNS = U 23 U 13 U 12 , we get 

/l \ /l 0\ 



u, 



23 



J- L 



, u 



13 



1 
\0 1/ 



(21) 



V° 71 75/ 

After rotating in the 23 and 13 plane, we can block-diagonalize the Mf 1 as 

/a + 2d V2b \ 



V2b a + b-d 
\ -a + b + dj 



U 12 = M, diag . 



(22) 



We can read m\ = \ — a+b+d\ 2 and obtain the remaining mass-squared eigenvalues, to 2 , m 2 ,, 
and the remaining mixing angle, 12 , by squaring (22), i.e. U^Nlf MfU MNS = (M^) 2 . 
There can be spontaneous CP violation and the parameters a, b, d are in general complex 
numbers: a = \a\e tiPa ,b = \b\e lVb , and d = \d\e tlfd . However, one of the vev's can be made 
real by SU(2) rotation. So we can set (p d = without loss of generality. Then it is 
straightforward to get the mixing angle 9 12 of Uu and the mass-squared eigenvalues: 



\d 
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™>i,2 = ^-|A=F A|, 
m 3 = \d\ 2 (x 2 + y 2 + 1 — 2xy cos ip a b — 2x cos ip a + 2y cos (fb) , 



A + E 

L o = tanfl 12 = \l - _ - , (23) 

where x = |a|/|d|, y = y? a 6 = y?a - y?6 and 

A = 2x 2 + 5y 2 + 5 + 2xy cos (fab + 2x cos (f a — 2y cos (fb, 
£ = — y 2 + 3 — 2xy cos Lp a b + 6rrcosy? a + 2y cos (fb, 



A = y^ + S^ + e 2 .), (24) 

where ei = y 2 + 2a;y cos ip a b + y cos and e 2 = 3y sin ipb- 

We can show that A > A (the equality sign holds for x = 27/32, y= (91- V4641) /128, 
and = (fib — ). So the smallness of the solar mass difference, Am| = Ara^ = |rf| 2 A, 
implies A rs 0. And this leads to E ps ei ~ e 2 ~ 0. From this we get a "magic relation" for 
y = 0(l): 

ipb ~ 0, y + 2a; cos y? a + 1 ~ 0, 
or ps 7r, y — 2x cos (p a — 1 rs 0. (25) 

We note that the above conditions give A w -3E or tan ^i 2 = l/>/2. Then we obtain the 
tribimaximal mixing (1) without further conditions 3 . Although the magic relation (25) 
is not guaranteed by any symmetry, it is satisfied in a large parameter space, especially 
when the CP violating phase ip a is allowed. In this approximation we get the following 



3 If (25) holds exactly, t\2 becomes undetermined. 
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mass-squared differences: 

Am^ = ml - m\ w \d\ 2 A, 

Am 2 ^ = ml — ml ps 6|d| 2 y (for = 0), or — 6|d| 2 ?/ (for = n). (26) 

Here Am^ > by definition and it can be identified with the Am 2 ol . However, the 
Am| 2 > can be either positive or negative: Am| 2 = Am 2 tm > for ip b ps (normal 
hierarchy) and Am| x = — Am 2 tm < for ipi, m n (inverted hierarchy). 

Now we get the numerical constraints on the five parameters, \d\,x,y,<p a and </?&, from 
the experimental data in Table 1. In Fig. 1, we show a scattered plot in the (ip a , plane. 
The blue (orange) color represents the case for the normal (inverted) hierarchy. We can 
see the allowed range of ipb is quite restricted and looks almost like line (its thickness is 
about 0.02). This shows the magic relation (25) works quite well. For the normal hierarchy 
((f b w 0) (f a is restricted in the region (tt/2,3tc/2), since the cosy9 a ps — (y + \)/2x < 0. 
by the magic relation. For the inverted hierarchy (<fb ps tt), since cosip a ~ (y — l)/2x, 
both signs are allowed in principle. However, the region in which both x and y are small is 
excluded and the minimum value of cos<p a allowed by the data is about —0.22. 

The very restricted range of fa signifies a fine-tuning. To avoid the fine-tuning problem, 
We set fa = 0, 71 identically for the normal and inverted hierarchy, respectively. To simplify 
expressions, we allow y to take negative values from now on, understanding that y > 
(y < 0) implies normal (inverted) hierarchy. Then the expressions for the mixing angle 612 
simplify greatly and give 

tan 2#i2 = ±^S|! for sign(y + 2x cos y? a + 1) = ±. (27) 

y ~ 3 

It is interesting to note that this expression is a function of y only and independent of x 
and (p a . The 3-a allowed range of y can be read from Fig. 2 and is given by 

y > 12.7 for (y + 2x cos ip a + 1) > 

1.14 < y < 1.70 for (y + 2x cos y? a + 1) < (28) 

for normal hierarchy and 

y < -4.71 for (y + 2xcos(p a + 1) < (29) 
9 



2.5 : J | : 

2.0 ; 1 1 : 

-Q " 1 1 

1.5 : i i : 

1.0; I I ; 

0.5 7 | l ; 

0.0 : 

1 2 3 4 5 6 

Figure 1: Scatted plot in the (ip a ,ip b ) plane. The blue (orange) color represents the case 
for the normal (inverted) hierarchy. The dashed black lines are ip a = 7r/2,37r/2. 

for inverted hierarchy. It is impossible to satisfy the data for the solar mixing angle when 
(y + 2x cos ip a + 1) > 0. 




Figure 2: Thick solid (dashed) curve: tan # 12 as a function of y for y+2x cos f a +l > (< 0). 
Dashed horizontal lines: 3-cr allowed range of tan# 12 , 0.577 < tan# 12 < 0.766. The left 
(right) panel is for normal (inverted) hierarchy. 

To consider the constraints from the mass-squared differences, we take the ratio 
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\Am 2 2 \/ because \d\ 2 is canceled in this case. Then Eqs. (23) and (24) give 

Arnjj \y + 2x cos ip a + 1| 



|Am| 2 | 



f(y) 



where 



f(y) 



^9-6?/ + V' 



(30) 



(31) 



Fig. 3 shows the plots of /(j/). 



2.0 



1.5 



0.5 - 



0.0 




10 



15 



\y\ 



Figure 3: Solid (Dashed) curve: f(y) for normal (inverted) hierarchy. 

The experimental value of p = (3.19 ± 0.19) x 10~ 2 can be satisfied by tuning y + 
2x cos ip a + 1 ~ 0. 

Using the relation, 2x cos ip a + y + 1 ~ 0, we can get 



_ _ y + i 

Am 2 tm ~~ G\y\ I V2cosy? a 



+ 2{y 2 -y+1) 



(32) 



in the leading order in p. The minimum of mi is obtained when cos 2 ip a = 1 for both 
normal and inverted hierarchy. Fig. 4 shows (mi) m 



Aml tm , that is, the mJWAm 2 



atni 



as a function of y when cos 2 ip a = 1. It shows that there is a lower bound on the lightest 
neutrino mass, muty > 0.03(0.05) eV for the normal (inverted) hierarchy. 
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-10 
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Figure 4: The mass ratio, (mi) min / a/ Am^ tm , as a function of y when cos 2 ip a = 1. The 
thick solid curves are allowed by the current experimental data. 

Now let us consider the bound on the effective Majorana mass for neutrinoless double- 
beta decay [4]. The amplitude of the neutrinoless double beta decay (Ou(3(3) is proportional 
to the effective Majorana mass for the Oz//5/5 defined by 

a + 2d\ = d\xe ilpa + 2| . (33) 

2(2/ (34) 

in the leading approximation in p. The ratio | (m^)\/ \J AmS| tm has lower bounds for the 
allowed range in y (thick parts in the solid blue curve in Fig. 5). Numerically we get 

\{m pp )\/y/Aml tm > 0.2 (for 1.13 < y < 1.7), > 0.57 (for y > 13), 
\(m^)\/VAml tm > 0.72 (for y < -4.5). (35) 

The model with A A symmetry which has non- vanishing lower bound for \(mpp)\/ a/ Am 2 tm 
even in the case of normal hierarchy was considered in [4]. Numerically we get similar 
values to theirs. 
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\( m f>p)\ 



m i (^MNs)li 



|(Mf) 



ii i 



Similar to (32), we get 



1(^)1 2 „ _L_ / ( y + i 

A ™atm %l V \2cOS(p a 




Am^ tm , as a function of y when cos 2 v? a = 1. The 



Figure 5: The mass ratio, (mpp) mm/ v ^.,,„ atm , 
blue (red, orange) curve corresponds to the case b = c (£ = 0.3, £ = —0.3). The thick parts 
are allowed by the current experimental data. 

3.2 The case b ^ c 

As mentioned in Section 2, the case b = c is not guaranteed by any symmetry. In this 
section we extend to the case b ^ c. The effective neutrino mass matrix given in (15) is 
now in the form: 

' a + 2d c b \ 

c b a — d 

b a — d c 



Ml 



off 



U{u)*M u U{uj 



(36) 



/ 



Since b = c case can already give the tribimaximal mixing matrix, we can see b and c 
cannot be so different. Therefore, we can apply the time-independent perturbation theory 
to diagonalize the neutrino mass matrix and expand in powers of £ = (c — b)/2b. Since b is 
assumed to be real, we also assume £ to be real. 

Since the 1st and 2nd eigenvalues giving the Am^ are quasi-degenerate, the blind 
application of the perturbation formula gives unreasonable results. To evade this problem 
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we diagonalized the 2x2 sub-matrix exactly. And then we applied the perturbation formula 
in the basis where the first two mass-squared eigenvalues are diagonal. 

The solar mass-squared difference is obtained to be 



Am^ ol ~ d 2 \2xcos(p a + y+l + £y|v / V - % + 9. (37) 

This implies the "magic relation" corresponding to (25) is simply replaced by 

2xcos(p a + y + 1 + £y « 0. (38) 

The atmospheric mass-squared difference is 

Aml tm ^6d 2 \y\(l + 0. (39) 

Since ( C 1, j/ > (?/ < 0) still gives normal (inverted) hierarchy. 

The correction in £ to the effective mass for the Of ft ft, (34), is given by 

K"W| 2 \ ( (Jl±±-) 2 -2(y-l) + 2yj V + 1 ) (40) 

The minimum values in tp a are obtained for <p a = 7r(0) for normal (inverted) hierarchy. 
These are plotted in Fig. 5 as a function of y for different values of £ = 0, ±0.3. The 
allowed regions are drawn in thick lines. The value |£| ~ 0.3 is almost maximum allowed 
by the 3-cr range in s 13 (see Figs. 6,7). Numerically the minimum values for the normal 
hierarchy are given by 



atm 



= 0.1, 0.2, 0.34 (for £ = 0.3, 0, -0.3, resp.). (41) 

a 

Since p <C 0.1, the correction to (40) in p does not change the results much. 
For the mixing angle 12 we get a formula similar to (27), 

2\/2\v\ 

tan 26»i2 ^ ± — for sign(y + 2x cos y? a + 1 + £y) = ±. (42) 

y-3 

The most significant change from the case b = c is that the non-vanishing #13 and 
consequently 5 for non-trivial ip a is allowed in b 7^ c case. The expression for s\ 3 = sin #13 
is obtained by 



s^^Jl + ( iy + 1 + i y)t ^ ) 2 . (43) 



V2 
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The CP violating phase S is given by 

r , -l f(y+ i + tan (p a 
6 ~ tan ' 



(mod 7r). 



(44) 



We do not have a definite prediction for s 13 . But the 3-cr range in Table 1 can be accom- 
modated. Fig. 6 (7) shows contours for the constant S13 and 5 for the normal (inverted) 
hierarchy case. We can see the 3-cr range for the S13 can be accommodated in the pertur- 
bative region for £. All the possible values of S are allowed by the current experimental 
values, although relatively small values of 5 are preferred. 

A very interesting prediction is for 623 ■ Up to the first order in £, the prediction for 623 
is still 7r/4. The first correction appears in second order in £: 



sin 2 623 = 0.5 + 0(£ 2 ) « 0.5 ± 0.01, 



(45) 



for the l-o" allowed £ (£ < 0.1) by S13. Therefore if the experiments would confirm significant 
deviation from 7r/4 for ^23, our model would be ruled out. 




3 Jj_p 0.0 3i 




2.0 2.5 3.0 3.5 4.0 4.5 



(fa 



(fa 



Figure 6: Solid (dashed) curves: contours for the constant S13 (5 (in degrees) ) for the 
normal hierarchy. The inner (outer) solid lines are contours for the 3-(1-)<t values of S13. 
We take y — 13 (y — 1.5) for the left (right) panel. 
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-0.4 
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Figure 7: The same with Fig. 6. But we take y = —5. 

4 Conclusions 

We studied a triplet Higgs model to generate Majorana neutrino masses and the mixing 
matrix in the framework of symmetry. With the assignments of A4 representations 
given in Table 2, we see that 

• The tribimaximal form of the neutrino mixing matrix can be naturally obtained for 
b = c. 

• There is a lower bound on the lightest neutrino mass: m\ > 0.03(0.05) eV for the 
normal (inverted) hierarchy. 

• There is a lower bound on the effective mass for the neutrinoless double beta decay: 
\{ m pp)\l a/ Am^tm > 0.2(0.57) for the normal (inverted) hierarchy. 

• For b 7^ c case, we can accommodate the data for the #13. 

• Even for b 7^ c case, the prediction for the atmospheric mixing angle does not change 
much from 6*23 = 7r/4 and gives sin 2 623 ~ 0.5 ± 0.01, which can be tested in near 
future. 
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